MASTER  COPY 


ke'ep  FOR  REPRODUCTION  PURPOSES 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  NO.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  easting  data  sources, 
gathenng  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comment  regarding  this  burden  estimates  cr  any  other  aspect  of  this 
collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  information  Operations  and  Reports,  1215  Jefferson 
Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20503. 


1 .  AGENCY  USE  ONLY  (Leave  blank) 


4.  TITLE  AND  SUBTITLE 


2.  REPORT  DATE 


3.  REPORT  TYPE  AND  DATES  COVERED 


Crystal  Chemistry  of  Ceramic/Metal  Systems 


6.  AU  i  fiCR(S) 

Robert  R.  Reeber 


-29  Sep  92  -  13  Dec  96 
5.  FUNDING  NUMBERS 


DAAL03-89-D-0001 


7.  PERFORMING  ORGANIZATION  NAMES(S)  AND  ADDRESS(ES) 

Department  of  Geology 

University  of  North  Carolina  at  Chapel  Hill 

Chapel  Hill,  NC  27599-33 15 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

U.S.  Armv  Research  Office 
P.O.Box '12211 

Research  Triangle  Park,,  NC  27709-2211 


10.  SPONSORING  /MONITORING 
AGENCY  REPORT  NUMBER 


AR0  30952.17-MS 


11.  SUPPLEMENTARY  NOTES 


The  views,  opinion*;  and/or  finding*;  contained  in  this  report  are  those  of  the  author(s)  and  should  not  be  construed  as 
an  official  Department  of  the  Army  position,  policy  or  decision,  unless  so  designated  by  other  documentation. 


12a.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 


12  b.  DISTRIBUTION  CODE 


Approved  for  public  release;  distribution  unlimited. 


13.  ABSTRACT  (Maximum  200  words) 

Robust  models  for  predicting  high  temperature/pressure  thermal  expansion  for  a  range  of  solids  have  been 
developed.  The  inter-relationships  between  several  thermophysical  properties  including  thermal  expansion,  specific  heat, 
elastic  properties,  and  self-diffusion  of  materials  has  been  investigated.  The  role  of  defects  on  thermal  expansion  at  high 
temperature  and  pressure  has  also  been  explored.  During  the  course  of  the  project  several  models  for  thermal  expansion  were 
developed  extending  earlier  work  of  Reeber.  The  product  of  thermal  expansion,  bulk  modulus,  and  volume  has  been  shown  to 
approach  a  constant  at  high  temperatures  for  alkali  halides,  rare  earth  oxides,  fee  metals,  and  ZnS.  Models  have  been  applied 
to  investigate 'high  temperature  thermal  expansion  of  alkali  halides,  rare  earth  oxides,  group  IV  elements,  fee  metals,  bcc 
metals,  SiC,  III-V  and  II- VI  compounds.  Analytical  expressions  have  been  provided  that  represent  the  thermal  expansion  of 
these  materials  from  near  0  K  to  the  melting  point.  Models  have  also  been  extended  to  high  pressure  thermal  expansion  of 
MgO,  NaCl,  Group  IV  elements,  and  W.  The  bulk  modulus  and  self-diffusion  for  W  and  Ta  have  been  predicted  to  their 
melting  points. 


19970210  198 


14.  SUBJECT  TERMS 


Thermal  Expansion,  High  Pressure,  High  Temperature,  Thermal  Defects 


15.  NUMBER  IF  PAGES 


16.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION 
OR  REPORT 

UNCLASSIFIED 


NSN  7540-01-280-5500 


18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

UNCLASSIFIED 


19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

UNCLASSIFIED 


20.  LIMITATION  OF  ABSTRACT 


Enclosure  1  DTjLC  i~vo 


Standard  Form  298  (Rev.  2-89) 

v  Prescribed  by  ANSI  Std.  239-18 
298-102 


Crystal  Chemistry  of  Ceramic/Mineral  Systems 


Final  Report 
Robert  R.  Reeber 

January  16, 1997 
U.S.  Army  Research  Office 
Contract  No:  DAAL03-89-D-0001 

30952-MS 

Department  of  Geology 
University  of  North  Carolina-CH 
Chapel  Hill,  NC  27599-3315 

APPROVED  FOR  PUBLIC  RELEASE; 

DISTRIBUTION  UNLIMITED. 


THE  VIEWS,  OPINIONS,  AND/OR  FINDINGS  CONTAINED  IN  THIS  REPORT  ARE 
THOSE  OF  THE  AUTHOR  AND  SHOULD  NOT  BE  CONSTRUED  AS  AN  OFFICIAL 
DEPARTMENT  OF  THE  ARMY  POSITION,  POLICY,  OR  DECISION,  UNLESS  SO 
DESIGNATED  BY  OTHER  DOCUMENTATION. 


A.  Research  Objectives: 


The  objectives  of  this  work  were  twofold  (1)  to  develop  x-ray  diffraction  techniques,  from 
an  experimental  and  analytical  standpoint,  for  Laue  and  continuous  x-ray  methods;  and  2)  to 
model  and  predict  physical  and  crystallochemical  properties  of  isostructural  inorganic  compounds 
with  emphasis  on  thermal  expansion,  phase  transitions  and  mechanical  properties. 

B.  Summary  of  Most  Important  Results: 

Significant  progress  was  made  on  analytical  interpretation  of  Laue  photographs  and  a 
Laue  program  written  in  Fortran  was  completed.  CdZnCaBaCu  perovskites  were  synthesized  by 
reacting  and  oxidizing  selenides  and  sulfides  of  Cd  and  Zn  with  a  CaBaCu  precursor.  The 
materials  were  evaluated  for  superconducting  properties  and,  although  structurally  related  to  high 
Tc  superconductors,  showed  no  indication  of  superconductivity  above  77  K.  A  single  phase 
CaBaCu  oxide  material  was  identified  as  orthorhombic  PMNA  space  group. 

The  main  emphasis  of  the  work  concentrated  on  Objective  2),  which  was  in  concert  with 
Mr.  Wang’s  doctoral  thesis.  This  objective,  to  model  and  predict  physical  and  crystallochemical 
properties  of  isostructural  inorganic  compounds,  was  so  successful  that  the  scope  was  extended 
to  include  refractory  metals  and  group  IV  semiconductors  including  diamond. 

Robust  models  for  predicting  high  temperature/pressure  thermal  expansion  for  a  range  of 
solids  have  been  developed.  The  inter-relationships  between  several  thermophysical  properties 
including  thermal  expansion,  specific  heat,  elastic  properties,  and  self-diffusion  of  materials  was 
investigated.  The  role  of  defects  on  thermal  expansion  at  high  temperature  and  pressure  was 
modeled  theoretically.  Exact  expressions  for  these  properties  were  provided  that  are  in  excellent 
agreement  with  existing  experimental  data. 

During  the  course  of  the  project  several  models  for  thermal  expansion  were  developed 
that  extended  earlier  work  of  Reeber.  The  product  of  thermal  expansion,  bulk  modulus,  and 
volume  has  been  shown  to  approach  a  constant  at  high  temperatures  for  alkali  halides,  rare  earth 
oxides,  fee  metals,  and  ZnS.  The  work  has  been  applied  to  model  and  predict  high  temperature 
thermal  expansion  of  alkali  halides,  rare  earth  oxides,  group  IV  elements,  fee  metals,  bcc  metals, 
SiC,  III-V  and  II-VI  compounds. 

Analytical  expressions  have  been  provided  that  represent  the  thermal  expansion  of  these 
materials  from  near  0  K  to  the  melting  point.  The  models  have  been  extended  to  predict  high 
pressure  thermal  expansion  of  MgO,  NaCl,  Group  IV  elements,  and  W.  The  bulk  modulus  and 
self-diffusion  for  W  and  Ta  have  been  predicted  to  their  melting  points.  Variation  of  thermal 
expansion  and  strength  with  pressure/temperature  for  diamond  and  tungsten  has  also  been 
studied.  One  report,  a  PhD  dissertation,  and  fourteen  peer  reviewed  papers  were  accepted  and 
published  during  the  course  of  the  research.  The  models  developed  will  have  impacts  on  both 
materials  science  and  mineral  physics. 
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Work  on  tungsten  and  diamond  was  presented  at  the  20th  Army  Science  Conference  and 
was  awarded  a  certificate  as  the  best  poster  paper  of  its  session.  The  bulk  moduli  predictions  for 
tungsten  and  tantalum  was  given  as  a  poster  paper  at  the  Army/University  of  Michigan 
Workshop  on  Hot  Gas  Erosion  of  Materials.  Requests  for  the  results  of  both  papers  were 
received  from  Dr.  John  Vasilakis,  Army  Large  Caliber  Research  Facility  at  Watervliet  Arsenal , 
Watervliet,  NY. 
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5.  Report  of  Inventions:  None 

Appendix  1  Abstracts  of  Papers  Listed  in  Part  C. 

1 .  The  role  of  surface  engineering  for  controlling  structural  ceramic  properties  is  reviewed. 
Initially,  a  perspective  of  the  effects  of  an  unmodified  surface  state  on  fracture,  adhesion, 
tribological  properties,  and  degradation  is  given.  Ceramic  surface  modification  techniques 

'  including  laser,  ion,  ion-beam-assisted  deposition,  and  plasma-assisted  deposition  are  then 
discussed  with  regard  to  how  each  non-equilibrium  technique  affects  bulk  and  surface 
properties.  The  importance  of  experiments,  process  and  materials  characterization,  and 
modeling  for  understanding  changes  in  localized  bonding,  chemistry,  and  structure  is 
emphasized.  Although  not  exhaustive,  a  broad  guide  to  important  literature  is  given. 

2.  Thermal  expansion  of  ionic  crystals  at  high  temperature  has  been  treated  as  contributions 
from  the  perfect  crystal  and  its  lattice  thermal  defects.  In  this  paper,  a  semi-empirical 
quasi-harmonic  model[l,2]  represents  the  thermal  expansion  of  the  perfect  crystal,  and 
contributions  from  lattice  defects  at  high  temperature  are  described  by  the  theory  of 
thermodynamics  of  point  defects.  The  formation  energies  of  thermal  defects  of  those 
crystals  considered  have  been  obtained  from  fitting  thermal  expansion  differences  between 
real  and  perfect  crystals.  Results  show  that  there  is  a  linear  trend  between  the  formation 
energies  of  thermal  lattice  defects  and  their  Debye  temperatures.  The  formation  energies 
obtained  are  smaller  than  theoretical  calculations  and  ionic  conductivity  measurements. 
This  indicates  thermal  defects  are  ordered  at  high  temperatures. 

3.  Low  temperature  caesium  halides  thermal  expansion  has  been  approximated  by  a  quasi- 
harmonic  approach  (QHM)  [Reeber(1975),  phys.  stat.  sol.  (a),  32,  321-331,  Reeber  & 
Haas(1984),  Thermal  Expansion,  8,  31-48]  and  predicted  at  high  temperature  with  a 
modified  quasi-harmonic  model  (MQHM)  [Wang  &  Reeber( 1994a),  submitted  to  J.  Phys. 
Chem.  Solids].  Formation  enthalpies  of  high  temperature  defects  have  also  been  obtained 
by  fitting  lattice  differences  between  a  real  and  an  approximate  perfect  crystal  (QHM). 
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Typical  values  of  the  coefficients  of  thermal  expansion  and  lattice  parameters  are  given. 

4.  Thermal  expansion  data  for  eleven  sodium  chloride  structure-type  solids  have  been 
evaluated  with  an  empirical  modified  quasi-harmonic  model  for  thermal  expansion.  Earlier 
it  was  shown  for  Group  IV  elements  and  10  Grimm-Sommerfeld  compounds  that  an 
empirical  quasi-harmonic  model  fits  experimental  data  from  nearly  zero  to  intermediate 
temperatures  in  the  vicinity  of  TL,  the  approximate  temperature,  where  all  lattice  modes 
are  fully  excited.  At  higher  temperatures,  experimentally  measured  thermal  expansion  is 
slightly  higher  than  an  extrapolation  of  this  empirical  relationship.  The  difference  can  be 
approximated  with  a  second-order  term  that  gives  excellent  quantitative  agreement  with 
experiment  for  the  available  data  to  approximately  90%  of  the  melting  temperature. 

5.  Thermal  expansion  and  molar  volume  of  periclase  are  important  parameters  for 
thermodynamic  modeling  of  the  earth's  interior.  We  have  critically  reviewed  the  available 
data  and  have  found  wide  variations  in  the  literature  for  MgO  lattice  parameters  at  room 
and  high  temperatures.  Earlier  Reeber  (1975)  developed  an  empirical  quasi-harmonic 
model  (QHM)  that  provides  an  excellent  representation  of  thermal  expansion  from  near  0 
to  1200K.  Wang  and  Reeber  (1995a)  show  for  alkali  halides  that  above  a  temperature  TL, 
the  approximate  temperature  where  all  lattice  modes  are  effectively  excited, 
experimentally  measured  thermal  expansion  is  slightly  higher  than  an  extrapolation  of  this 

,  empirical  relationship.  This  difference  is  approximated  with  a  2nd  order  term  that  gives 
excellent  quantitative  agreement  with  experiment  for  alkali  halides  and  cesium  halides 
(Wang  and  Reeber,  1995b)  to  approximately  90%  of  the  melting  temperature.  The 
existing  literature  for  periclase  is  evaluated  within  the  constraints  of  the  modified  quasi¬ 
harmonic  model  (MQEOVI)  and  x-ray  measurements  at  room  temperature.  At  very  high 
temperatures  where  experimental  data  is  not  available  the  MQHM  provides  a  method  to 
estimate  High  temperature  thermal  expansion.  Our  best  estimate  of  lattice  parameter  at  OK 
is  4.20565 A  and  at  298K  4.2 1 1 54A. 

6.  Thermal  expansion  of  MgO,  periclase,  is  approximated  by  contributions  from  a  quasi¬ 
harmonic  perfect  crystal  and  its  equilibrium  high  temperature/pressure  defects.  We 
demonstrate  that  the  influence  of  thermal  defects  on  high  pressure  thermal  expansion  is 
very  small  and  may  be  neglected  at  lower  mantle  conditions.  Following  that,  a  simple 
relationship  has  been  obtained  between  high  pressure  thermal  expansion  and  the  thermal 
expansion  of  a  quasi-harmonic  perfect  crystal  at  zero  pressure.  Our  estimation  of  MgO 
thermal  expansion  at  the  core-mantle  boundary  is  9.36xlO'6K‘1. 

7.  Thermal  expansion  is  important  for  predicting  residual  stresses  in  epitaxial  films, 
composites  and  electronic  devices  as  well  as  for  providing  information  relevant  to  an 
understanding  of  interatomic  potentials  and  the  equation  of  state  of  materials.  Model 
calculations  have  many  assumptions,  both  inherent  and  implicit,  and  have  difficulty 
accurately  representing  thermal  expansion  at  high  temperatures  and  pressures.  We  utilize  a 
semi-empirical  quasi-harmonic  model  to  evaluate  available  data  for  P-silicon  carbide, 
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gallium  phosphide  and  indium  phosphide.  The  model  allows  prediction  of  the  thermal 
properties  of  these  semiconductors  from  near  o  K  to  the  vicinity  of  their  melting  points. 
The  approach,  consisting  of  a  simplified  frequency  spectrum  with  several  Einstein  terms, 
provides  a  convenient  mathematical  method  where  a  minimum  of  empirical  parameters 
represent  the  thermal  property. 

8.  Thermal  expansion  is  practically  and  theoretically  important  for  predicting  residual  stresses 
in  epitaxial  films  and  electronic  devices  as  well  as  for  providing  information  about  the 
anharmonic  properties  of  materials,  and  determining  their  equation  of  state.  Model 
calculation  have  had  difficulty  in  accurately  representing  thermal  expansion  at  high 
temperatures  where  such  considerations  are  most  important.  We  utilize  a  semi-empirical 
quasi-harmonic  model  to  evaluate  available  data  for  silicon  and  germanium.  The  model 
allows  us  to  predict  the  thermal  properties  of  these  semiconductors  from  near  OK  to  the 
vicinity  of  their  melting  points.  The  approach,  consisting  of  a  simplified  frequency 
spectrum  with  several  Einstein  terms,  provides  a  convenient  mathematical  method  where  a 
minimum  of  empirical  parameters  represent  the  thermal  property. 

9.  Thermal  expansion  is  of  both  practical  and  theoretical  importance  and  it  together  with 
specific  heat  is  essential  for  predicting  a  thermodynamic  equation  of  state.  We  utilize  a 
semi-empirical  quasi-harmonic  model  to  evaluate  available  data  for  diamond.  The  model 

"  ,  allow  us  to  predict  the  thermal  properties  of  the  metastable  diamond  polymorph  to  3000K. 
The  approach  consisting  of  a  simplified  frequency  spectrum  with  several  Einstein  terms 
provides  a  convenient  mathematical  method  where  a  minimum  of  empirical  parameters 
represent  the  thermal  property. 

10.  The  product  of  thermal  expansion,  bulk  modulus,  and  volume,  a^rF,  is  the  partial 
temperature  derivative  of  the  work  done  by  thermal  pressure.  For  copper  this 
thermodynamic  product  resembles  a  specific  heat  (Cv)  curve  and  approaches  a  constant  at 
high  temperature.  Our  recent  developed  model  utilizes  available  data  for  this  parameter 
and  recommends  new  copper  thermal  expansion  values  from  0  K  to  the  melting  point. 

1 1 .  Thermal  defects  in  ionic  materials  can  have  important  effects  on  their  thermal  expansion  at 
high  temperatures.  Earlier  treatments  of  thermal  expansion  generally  have  neglected  or  not 
considered  such  effects.  Here  an  analytical  expression  for  the  thermal  defect  contribution 
ancl  its  dependence  on  pressure  is  derived.  We  show  that  such  contributions,  whcih  are 
significant  at  high  temperature  and  atmospheric  pressure,  become  negligible  at  pressures 
above  approximately  0.25  to  0.35  of  the  bulk  modulus  at  standard  conditions.  At  very 
high  pressure,  based  on  Birch's  (1968)  relationship  between  high  and  low  pressure  thermal 
expansion,  and  assuming  a KT  is  independent  of  pressure,  NaCl  thermal  expansion  can  be 
calculated  within  the  constraints  of  a  semi-empirical  quasi-harmonic  perfect  crystal  model. 
The  calculations  are  compared  with  available  theoretical  and  experimental  values  over  an 
extented  temperature/pressure  regime.  The  method  should  have  broad  applicability  for 
other  ionic  crystals. 
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12.  In  this  paper,  a  new  set  of  experimental  data,  avKTV,  representing  the  partial  temperature 
derivative  of  the  work  done  by  the  thermal  pressure  of  solid,  is  fitted  by  n-terms  of  a 
modified  Einstein  model.  Experimental  data  show  that  CLyKTV,  not  avKT,  approaches  a 
constant  value  at  high  temperature.  Based  on  the  observed  linear  relationship  of  isothermal 
bulk  modulus  with  temperature  at  high  temperature,  thermal  expansion  can  be  evaluated 
by  fitting  apXTV  data.  Our  previous  results  have  shown  that  at  low  temperature  or  for 
materials  with  less  variable  bulk  modulus  and  expansivity,  thermal  expansion  data  can  be 
simply  approximated  by  an  n  term  Einstein  model.  More  generally  and  for  many  materials, 
avKTV  data  resembles  an  isochoric  specific  heat  curve.  With  this  method,  thermal 
expansion  can  be  predicted  at  high  temperature  from  low  and  intermediate  temperature 
range  data.  With  accurate  thermal  expansion  data,  high  temperature  bulk  moduli  can  also 
be  predicted. 

1 3 .  Thermal  expansion  at  high  temperatures  and  pressures  is  an  excellent  predictor  of  ultimate 
strength.  Theoretical  treatments  have  not  been  successful  for  accurately  predicting  thermal 
expansion  at  high  temperatures.  Here  we  approximate  the  thermal  expansion  of  group  IV 
semiconductors  and  tungsten  by  contributions  from  a  quasi-harmonic  perfect  crystal  and 
its  equilibrium  high  temperature/pressure  defects.  At  higher  pressures  thermal  defect 
effects  are  small  and  may  be  neglected.  The  thermal  expansion  for  each  element  is 
calculated  from  OK  to  its  respective  melting  point  and  for  range  of  pressure  from  0  GPa  to 

'  ,s  1 50  GPa.  It  is  anticipated  that  an  understanding  of  thermal  defect  behavior  as  a  function  of 
temperature  and  pressure,  when  contrasted  with  equivalent  shock  experiments,  will 
provide  insights  for  modeling  and  design  of  improved  penetrator  materials. 

14.  The  tungsten  bulk  modulus  is  calculated  from  self-diffusion,  thermal  expansion  and 
specific  heat  data.  Our  estimate  of  the  adiabatic  bulk  modulus  as  a  function  of  temperature 
is  about  8%  higher  than  earlier  results  at  the  highest  temperature.  For  tantalum,  self¬ 
diffusion  coefficients  are  calculated  from  high  temperature  thermoelastic  data. 

Comparisons  with  available  experimental  results  and  earlier  calculations  are  provided. 
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